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Abstract 

We attempt to simultaneously explain the recently observed 3.55 keV X-ray line in the analysis 
of XMM-Newton telescope data and the galactic center gamma ray excess observed by the Fermi 
gamma ray space telescope within an abelian gauge extension of standard model. We consider 
a two component dark matter scenario with tree level mass difference 3.55 keV such that the 
heavier one can decay into the lighter one and a photon with energy 3.55 keV. The lighter dark 
matter candidate is protected from decaying into the standard model particles by a remnant Z -2 
symmetry into which the abelian gauge symmetry gets spontaneously broken. If the mass of the 
dark matter particle is chosen to be within 31 — 40 GeV, then this model can also explain the 
galactic center gamma ray excess if the dark matter annihilation into bb pairs has a cross section 
of (av) — (1.4 — 2.0) x 10 -26 cm 3 /s. We constrain the model from the requirement of producing 
correct dark matter relic density, 3.55 keV X-ray line flux and galactic center gamma ray excess. We 
also impose the bounds coming from dark matter direct detection experiments as well as collider 
limits on additional gauge boson mass and gauge coupling. We also briefly discuss how this model 
can give rise to sub-eV neutrino masses at tree level as well as one-loop level while keeping the dark 
matter mass at few tens of GeV. We also constrain the model parameters from the requirement 
of keeping the one-loop mass difference between two dark matter particles below a keV. We find 
that the constraints from light neutrino mass and keV mass splitting between two dark matter 
components show more preference for opposite CP eigenvalues of the two fermion singlet dark 
matter candidates in the model 
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I. INTRODUCTION 


Recent analysis [1, 2] of the observations made by XMM-Newton X-ray telescope have 
pointed towards a monochromatic X-ray line with approximate energy 3.55 keV in the spec¬ 
trum of 73 galaxy clusters (For a review of dark matter, please see [3]). The same line also 
appears in the Chandra observations of the Perseus cluster [1], In the absence of any astro- 
physical interpretation of the line due to some atomic transitions, the origin of this X-ray 
line can be explained naturally by sterile neutrino dark matter with mass approximately 
7.1 keV decaying into a photon and a standard model (SM) neutrino. This was pointed out 
by the authors in [1, 2] and subsequently studied within the framework of specific models 
[4]. Several other particle physics explanations of the X-ray line have also been put forward 
in [5, 6]. Although most of the particle physics explanations consider late time decay or 
annihilation of multi-keV dark matter particles as the origin of the X-ray line, there have 
also been a few discussions on the scenario where the X-ray line can be generated by transi¬ 
tions between electroweak scale dark matter states with keV mass splittings [7, 8]. In spite 
of the fact that both keV scale as well as weak scale dark matter candidates can explain 
the same signal, their implications in cosmology and astrophysical structure formation can 
be very different. As stated in [1, 2], a keV scale sterile neutrino should have mixing with 
the SM neutrinos of the order « 1CU 11 — 1CT 10 to explain the X-ray line. Such a tiny mix¬ 
ing prevents the sterile neutrinos from entering thermal equilibrium in the early Universe, 
making it necessary to have some additional physics responsible for the production of sterile 
neutrinos. However, electroweak scale dark matter particles can be thermally populated 
in the early Universe due to their sizable interactions either through gauge bosons, Higgs 
portals or fermion portals etc. These scenarios are studied in the context of so-called weakly 
interacting massive particle (WIMP). The fact that the WIMP annihilation cross section 
turns out to be almost equal to the annihilation cross section of thermal dark matter in order 
to produce the correct dark matter relic abundance observed by the Planck experiment [9] 

Q DM h 2 = 0.1187 ±0.0017 (1) 

is known as the WIMP Miracle. In the above equation (1), is the density parameter and 
h = (Hubble Parameter)/100 is a parameter of order unity. 

Motivated by the possibility of explaining the origin of 3.55 keV X-ray line within WIMP 
dark matter framework [7, 8], here we consider an abelian gauge extension of SM with two 
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Majorana fermion dark matter candidates with a keV mass splitting. The UV complete 
model, originally proposed by [10] and later studied in the context of dark matter and eV 
scale sterile neutrino in [11] and [12] respectively, can naturally explain dark matter and the 
origin of tiny neutrino masses. Sub-eV scale SM neutrino masses arise both at tree level as 
well as one-loop level with dark matter particles running inside the loops, a framework more 
popularly known as "scotogenic" model [13]. Recently, this model was also studied [14] in the 
context of explaining the galactic center gamma ray excess observed by the Fermi Gamma 
Ray Space Telescope [15]. The abelian gauge charges of the SM as well as beyond SM fields 
are chosen in such a way that the model is free from gauge anomalies and the abelian gauge 
symmetry U(l)x gets spontaneously broken down to a remnant Z 2 symmetry so that the 
lightest Z 2 -odd particle is stable and hence can be a dark matter candidate. As studied 
in details in [11], this model has several dark matter candidates namely, fermion singlet, 
fermion triplet, scalar singlet and scalar doublet. Scalar dark matter phenomenology is 
similar to the Higgs portal models discussed extensively in the literature. In these scenarios, 
the scalar dark matter self-annihilates into the Standard Model (SM) particles through the 
Higgs boson. Co-annihilations through gauge bosons can also play a role if the CP even and 
CP odd components of the neutral Higgs have a tiny mass difference as discussed recently 
in [16] within the context of a different model. 

Instead of pursuing Higgs portal like scalar dark matter scenarios in the model, we study 
the fermionic dark matter sector. This is also relevant to our discussion on the origin of 
3.55 keV X-ray line. This is because, if transition between two semi-degenerate weak scale 
dark matter candidates with keV mass splitting is the origin of the X-ray line, then the 
dark matter candidates have to be fermions as one scalar decaying into another scalar and 
a photon does not conserve spin. As we will see in the next section, our model has two 
fermion singlet dark matter candidates with different gauge charges and two fermion triplet 
dark matter candidates with the same gauge charge. We can choose either a triplet-singlet 
or a singlet-singlet combination of two semi-degenerate dark matter candidates to explain 
dark matter abundance as well as the origin of 3.55 keV X-ray line simultaneously. However, 
the neutral component of fermion triplet needs to be very heavy (2.28 — 2.42 TeV) in order 
to reproduce correct dark matter relic density [17]. To allow the possibility of low mass dark 
matter, we therefore confine our discussion to fermion singlet dark matter in this work. That 
is, we explore the possibility of two fermion singlet dark matter candidates with keV mass 


4 


splitting in this model which can simultaneously give rise to the 3.55 keV X-ray line and 
satisfy experimental bounds on dark matter relic density as well as direct detection cross 
section. Such fermion singlet dark matter particle will self-annihilate through the abelian 
vector boson X into SM particles. We also incorporate the collider constraints on such 
additional vector boson and its gauge coupling. We find that, although the relic density and 
direct detection constraints allow a significant region of the parameter space, the collider 
constraints reduce the parameter space into the s-wave resonance region where the gauge 
boson mass is approximately twice that of dark matter mass. We constrain the model further 
by incorporating the bound from X-ray line data on the decay width of the heavier dark 
matter particle. We then check whether the same model can also give rise to the galactic 
center (GC) gamma ray excess observed by the Fermi Gamma Ray Space Telescope which 
has a feature similar to annihilating dark matter [15]. Finally, we briefly discuss whether 
the chosen dark matter masses are compatible with sub-eV SM neutrino masses and also 
constrain the model parameters in order to keep the one-loop mass splitting between two 
dark matter particles below keV scale. 

This paper is organized as follows: in section II, we briefly discuss the model. In section 

III, we discuss two component dark matter scenario as a source of 3.55 keV X-Ray line and 
GC gamma ray excess taking into account all necessary experimental constraints. In section 

IV, we discuss the compatibility of light singlet fermion dark matter with neutrino mass and 
in section V, we discuss the one-loop mass splitting between two dark matter candidates. 
Finally, we conclude in section VI. 


II. THE MODEL 

The model which we take as a starting point of our discussion was first proposed in [10] 
which has fermion content shown in table II. The scalar content of the model is modified in 
our present work with a different U(l)x gauge charge for the singlet scalar \2 and a newly 
added scalar singlet X 5 as shown in table II. 

The third column in table II shows the U{l)x gauge charges of various fields which 
satisfy the gauge anomaly matching conditions. The charges of the scalar fields in table II 
are chosen according to the desired neutrino and dark matter phenomenology. The Higgs 
content chosen in the model is not arbitrary and is needed, which leads to the possibility of 
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Particle 

SU{ 3) c x SU(2) l x U{l) Y 

U( l)x 

^2 

(■ u,d) L 

(3,2, |) 

Tl\ 

+ 

UR 

(3,1,|) 

1 

CO 

3 

4 ^ 

+ 

dR 

(3,1,4) 

|(ni + 3n 4 ) 

+ 

{v, e) L 

(1.2,4) 

77-4 

+ 

cr 

(1,1,-1) 

l(-9ni + 5n 4 ) 

+ 

Nr 

( 1 , 1 , 0 ) 

§(3ni + n 4 ) 

- 

'ZlR,2R 

(1,3,0) 

|(3ni + n 4 ) 

- 

Sir 

( 1 , 1 , 0 ) 

l(3ni + n 4 ) 

+ 

S 2 R 

( 1 , 1 , 0 ) 

-|(3ni + n 4 ) 
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TABLE I: Fermion Content of the Model 


Particle 

SU( 3)c x SU(2) l x U(l) Y 

U(l)x 

^2 

(<£ + ,0°)l 

(1,2,-|) 

|(ni - n 4 ) 

+ 

W> + ,/)2 

(1,2,-i) 

l(9ni - n 4 ) 

+ 

(^ + ,/) 3 

(1,2,4) 

|(9ni - 5n 4 ) 

- 

Xi 

( 1 , 1 , 0 ) 

— ^(3ni + n 4 ) 

+ 

X2 

( 1 , 1 , 0 ) 

-|(3ni + n 4 ) 

+ 

X3 

(1,1,0) 

— |(3ni + n 4 ) 

- 

X4 

(1,1,0) 

— f(3ni + n 4 ) 

+ 

X5 

(1,1,0) 

(3n 4 — n 4 ) 

+ 


TABLE II: Scalar Content of the Model 


radiative neutrino masses as well as a remnant Z 2 symmetry in a manner proposed in [13]. 
In this model, the quarks couple to scalar $1 and charged leptons to the scalar $2 whereas 
(y,e)L couples to Nr,Zr through <h 3 and to Sir through dq. The Lagrangian which can 
be constructed from the above particle content has an automatic Z 2 symmetry and hence 
the model provides a stable cold dark matter candidate in terms of the lightest odd particle 
under this Z 2 symmetry. The Z 2 transformations of the fields are shown in the fourth column 
of table II and table II. The scalar Lagrangian relevant for future discussion can be written 
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as 


K D f3Xixi®\®3 + hxlx4^l^2 + f6{®\®l)X3X5 (2) 

Similarly, the relevant part of the Yukawa Lagrangian for the model can be written as 

Cy D yL&\SiR + hjyL&lN R + IiyL&1'E r + fNN R N R x 4 + fsS\ rS ]R xi 

+ f^R^RX4 + fS 2 S 2 RS 2 RXJ 2 + fl2SlRS 2 Rxl (3) 

Let us denote the vacuum expectation values (vev) of various neutral scalar fields as (0° 2 ) = 
v i, 2 , (x? 245 ) = 'tii, 2 , 4 , 5 ■ The additional U(l)x gauge boson mass, which is relevant for dark 
matter phenomenology is given by 

XM 2 , 1 

Mx = 2g 2 x (-—^(9n 1 -n /L )(n 1 -7i 4 ) + — (3n 1 +ru)' 2 (4u 2 1 + 25ul + 9ul) + (3ni-n4) 2 ul) (4) 
*9 2 -Lu 

where g\ is the U(l)x gauge coupling. For simplicity, the mixing between the neutral 
electroweak gauge bosons and the additional U(l)x gauge boson is chosen to be zero which 
gives rise to the following constraint 

3(n 4 — n i)vl = (9ni — n 4 )v \ (5) 

which further implies 1 < n 4 jn\ < 9. 


III. SINGLET FERMION DARK MATTER 

The relic abundance of a dark matter particle x i s given by the Boltzmann equation 

+ 3Hn x = -W v )( n l - K q ) 2 ) (6) 

where n x is the number density of the dark matter particle y and n® q is the number density 
when x was in thermal equilibrium. H is the Hubble expansion rate of the Universe and (av) 
is the thermally averaged annihilation cross section of the dark matter particle y. In terms 
of partial wave expansion (av) = a + bv 2 . Numerical solution of the Boltzmann equation 
above gives [18] 

2 1.04 x 10% F 

x M PI Jg;(a + 3b/x f ) [> 

where Xf = m x /Tf, Tf is the freeze-out temperature, g * is the number of relativistic degrees 
of freedom at the time of freeze-out. Dark matter particles with electroweak scale mass 
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and couplings freeze out at temperatures approximately in the range corresponding to Xf ~ 
20 — 30. More generally, xj can be calculated from the relation 


x f 


In 


0.038gMpim x < crv > 


l/2„ 


.1/2 
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( 8 ) 


where g is the number of internal degrees of freedom of the dark matter particle x and Mpi 
is the Planck mass. The thermal averaged annihilation cross section (crv) is given by [19] 


(av) = 


8m*TK$(m x /T) J Am a 


cr(s — 4m 2 x )y/sKi(y/s/T)ds 


(9) 


where K, 's are modihed Bessel functions of order i, m x is the mass of Dark Matter particle 
and T is the temperature. 

There are two singlet fermions Nr, S 2 R in this model which are odd under the remnant 
Z 2 symmetry and hence can be a dark matter candidate, if lightest among all the Z 2 -odd 
particles. We consider a scenario where S 2 R is the lightest and Nr is the next to lightest 
Z 2 -odd particle of the model. If the lifetime of Nr is very high, longer than the present 
age of the Universe, then both S 2 R and Nr can contribute to the present abundance of 
dark matter. From the held content and their gauge charges, one can see that there is no 
term in the Lagrangian which involves both S 2 R and Nr. Also there is no scalar which 
couples to both S 2 R and Nr. Thus, there is no co-annihilating processes between S 2 R and 
Nr which can contribute to the dark matter relic abundance. Hence, one can calculate the 
relic abundance of S 2 r and N R separately, keeping them decoupled. To calculate the relic 
density of either S 2 r or Nr, we need to find out its annihilation cross-section to standard 
model particles. For zero Z — X mixing, the dominant annihilation channel is the one with 
X boson mediation. Since the singlet fermions are of Majorana type, they have only axial 
coupling to the vector boson. The annihilation cross-section of Majorana singlet fermion 
into SM fermion anti-fermion pairs // through s-channel X boson [20] can be written as 

1 — 4 m 2 /s 1 1//2 


n. 


a = 


127rs [(s — m x ) 2 + M X T X ] 
6 s 


2 2 

9fa9 X a 


4 mi 


m 


7- 


M% 


+ 


1 - 4 M 2 x /s 
3s 2 


x 


M A X 


— s 


+ s(s — 4m 2 )^ 


+ 9%9 2 xa ( s + 2 m 2 f )(s - 4m 2 ) 


( 10 ) 














Expanding in powers of v 2 gives av in the form a + bv 2 where a and b are given by 


a = 


ncg 2 fa rn 2 f gl a m 2 x 


24? T 2 m 2 MM 2 x - Am 2 )' 2 + M 2 X Y 2 X ) 
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m f 
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mt 
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M\ 


+ 192 
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m| 
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2m 2 (M| — 4m 2 
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4 (M| — 4m 2 ) 2 + M 2 T 2 X 8(m 2 — m 2 )m 2 
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9 fa 
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ml 
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( 11 ) 


Lf ±V± X 1V1 X ^ 

The Decay width of the X boson denoted by Tx is given by 

3/2 


T*-xx 

r A'^/7 


n c M x g\ 

Yl-nS 


Am 2 
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n r M 


x 


12itS 


1 - 


Am 2 f 

mi 


1/2 


9% 


1 - 


4m^ 


m 


x 


+ g} v I 1 + 2 


mi 


'Ml, j 


( 12 ) 


The mass of the gauge boson X in the above expressions is given by equation (4). For 
simplicity, we assume u\ = U 2 = = u such that the mass of X boson can be written as 


M‘x — 2 g\ 


3 m w 


8 g. 


19 

~8~ 


2 (9n 4 — n 4 )(n 4 — n 4 ) + 4(3ni + nf) 2 u 2 + (3n 4 — n 4 ) 2 w 2 


(13) 


The couplings gf v ,gfa, 9xvi 9\a of fermions and dark matter to X boson are tabulated in the 
table III. 



n c 

9fv/gx 

9fa/gx 

l = e, fi,T 

1 

| (n 4 - m) 

| (n 4 — 9n 4 ) 

Vl 

1 

714 

2 

714 

2 

U = u,c 

3 

|(llni — n 4 ) 

e 

1 

coloo 

D = d,s,b 

3 

g(5n 4 + 3n 4 ) 

|(n 4 - n 4 ) 

Nr 

1 

0 

|(3n 4 + n 4 ) 

S 2 R 

1 

0 

-|(3ni + n 4 ) 


TABLE III: Couplings of SM particles and dark matter to the vector boson X 


Using the couplings given in table III, we first check whether there exists a freeze-out 
temperature Tf for singlet fermionic dark matter such that for T > Tf, the singlet fermions 
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nriDM = 35 GeV; M x = 70 GeV; g x = 0.005 



m DM /T 


FIG. 1: Comparison of self-interaction rates of dark matter candidates (Nr, S2r) with ms 2R = 
ttidm = 35 GeV and mx R ~ U7g 2i? = 3-55 keV with the Hubble expansion rate H for Mx = 70 GeV, 
gx = 0.005. 

enter thermal equilibrium. This is crucial in order to use the standard relic abundance 
formula given by equation (7) for WIMP dark matter. For this we need to calculate the 
annihilation rate of singlet dark matter particles and compare with the Hubble expansion 
rate of the Universe. To calculate the interaction rate and hence annihilation cross section 
in our case, we fix the gauge charges 771,774 but vary the U(l)x gauge coupling gx and 
gauge boson mass Mx- Similar to our approach in [11, 14], here also we choose a specific 
value of 77 ! from which 774 can be determined using the normalization relation n\ + n\ = 1. 
Using the same normalization, the 90% confidence level exclusion on Mx/gx was shown 
in [ 10 ] where the lowest allowed value of Mx/gx was found to be approximately 2 TeV 
for (f) = tan _ 1 ( 774 /? 7 i) = 1.5. Using this and the normalization relation involving 774,774 we 
determine both 774 , 774 . Since U(l)x gauge charges of all the fields are written in terms of 
774,774 it is sufficient to choose just these two values to determine all the gauge charges. After 
fixing dark matter mass as well as 774 4 , we vary gx and u and compute the annihilation 
cross section of dark matter particles. The freeze-out temperature Tf is then calculated 
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0.1 


<|>=tarr 1 n 4 /n 1 =1.45 



M x /g x > 2 TeV 


M x /g x > 3 TeV 


M x /g x >4 TeV 


Galactic Center Bound 


Relic Bound 

X 

XENON 100 bound 
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FIG. 2: Parameter space in the gx ~ Mx plane for two-component dark matter (Nr, S^b) scenario 
with mg 2R = 35 GeV and mx R ~ m S 2R = 3-55 keV. The red-hatched, green and blue dot-dashed 
regions correspond to the allowed region after the constraints on Mx/gx are imposed. The area 
to the left of the black line is ruled out by XENONIOO bounds on direct detection cross section. 
The solid red region corresponds to the parameter space favored by the relic density constraint. 
The solid black region corresponds to the parameter space favoured by galactic centre gamma ray 
excess. 

numerically by using the equation 

. 0.038 gm PL m Y < av > 

pi — m-—- 

c 111 1/2 1/2 
9* x f 

This is a simplihed form of equation (8). For a fixed value of dark matter mass m x , the 
annihilation cross section cr depends upon gx, Mx- For a particular pair of gx and Mx, we 
use this value of Xf and compute the relic abundance using equation (7). 

To show that the dark matter candidates in our model have gone through this generic 
freeze-out process, we compare their interaction rates with the Hubble expansion rate of 
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FIG. 3: Parameter space in the g\ — u plane for two-component dark matter (Nr, S^r) scenario 
with rns 2R = 35 GeV and mx R ~ m S 2 R = 3-55 keV. The red-hatched, green and blue dot-dashed 
regions correspond to the allowed region after the constraints on Mx/gx are imposed. The area 
to the left of the black line is ruled out by XENON100 bounds on direct detection cross section. 
The solid red region corresponds to the parameter space favored by the relic density constraint. 
The solid black region corresponds to the parameter space favoured by galactic centre gamma ray 
excess. 


the Universe. The interaction rate is given by T = n(av ), where n is the number density, 
and ( av ) can be calculated using equation (9). For a non-relativistic dark matter particle 
of mass m, the equilibrium number density is given by 

3/2 


n = g 


mT 

~2 7 


D -m/T 


(15) 


where we have taken the chemical potential to be zero. Here, g = 2 for Majorana fermion 
dark matter. On the other hand, the Hubble expansion parameter for the early radiation 
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FIG. 4: Relative contribution of Nr and S 2 R to dark matter relic abundance for rns 2R = 35 GeV 
and m Xr — ms 2 R = 3-55 keV. 


dominated epoch can be written as 


H = 1.66,!/^ («) 

N1 pi 

We plot T = n(av) as well as H as a function of m/T for two dark matter particles with 
masses rris 2R = m DM = 35 GeV and rri Xlt — m S2R = 3.55 keV for gauge boson mass M x = 70 
GeV and gauge coupling g x = 0.005. This is shown in figure 1. The point at which the 
interaction rate T falls below the Hubble expansion rate H corresponds to the freeze-out 
temperature Tf. From figure 1, we see that this crossover occurs at itidm/T ~ 15, which 
corresponds to freeze-out temperature Tf fa 2.33 GeV. Calculation of dark matter self- 
annihilation cross section also allows us to find out the parameter space giving rise to the 
correct relic abundance. This parameter space in g x — M x and g x — u planes are shown in 
2 and 3 respectively. 

We then consider the bound on dark matter nucleon scattering cross section from direct 
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detection experiments. Since both the dark matter candidates in our model are Majorana 
fermions, the vector current vanishes and they give rise to spin dependent scattering cross 
section with nuclei. However, spin independent scattering can arise if there exist scalar 
mediated interactions between dark matter and nucleons. In the limit where the mixing 
between singlet scalars X' 2 , X4 and the scalar doublet dq is negligible, we can consider only 
the gauge boson mediated spin dependent scattering between dark matter and nucleons. The 
latest upper bound on this scattering cross section comes from the XENONIOO experiment 
[21]. The expression for this spin dependent scattering of dark matter particles off nuclei 
through t-channel mediation of X boson can be written as 



where 


m x m N 

l^x N = j . t - 
m x + my 

and Jtv is the spin of the Xenon nucleus. The standard values of the nuclear quark content 
areA,^ = = 0.84 and A= A^ p) = —0.43 [22]. The average spins (S p ) and (S n ) of 

the Xenon nucleus are taken from [21] as given in table IV. The lowest upper bound at 90% 


Nucleus 

(Sn) 

(S P ) 

129 Xe 

0.329 

0.010 

131 Xe 

-0.272 

-0.009 


TABLE IV: Average Spin of Xenon Nucleus 


confidence level from XENONIOO experiment on spin dependent dark matter nuclei cross 
section is 3.5 x 10 -40 cm 2 for dark matter mass of 45 GeV. Here we take this as a conservative 
upper bound on direct detection cross section and show the region of parameter space in 
both gx — Mx and gx — u planes which gives rise to this cross section. This gives rise 
to a solid exclusion line in the figure 2 and 3 so that the region of parameter space above 
or towards left of this line is ruled out. Similar to the discussion in our earlier work [14], 
we also incorporate the collider bounds on Mx and gx- Collider constraints on additional 
gauge bosons masses with generic SM like gauge couplings force them to be heavier than 
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approximately 2.5 TeV. However, as discussed in [23], the bounds on the mass of additional 
boson X can be relaxed if it has non-negligible coupling to the dark sector. The authors 
showed that for X decaying into SM particles with branching ratio 90% and gx = 0.1, the 
lowest allowed value of Mx/gx is approximately 2.6 TeV. This limit goes up to 4 TeV and 4.4 
TeV, if gx is increased to 0.3 and weak gauge coupling g respectively. To apply these bounds, 
we calculate the branching ratios of X boson into SM and dark sector particles and find that 
the maximum branching ratio of X boson into dark matter particles is approximately 8.5%. 
According to the analysis of [23], this will correspond to an approximate bound Mx/gx > 2.6 
TeV for gx = 0.1. However, these bounds will be weaker if gx is lowered down into the 
resonance region that can be seen from figure 2 and 3. We apply moderate as well as 
conservative bounds on Mx/gx between 2 TeV to 4 TeV and show the portion of parameter 
space left after that in figure 2 and figure 3. 

It can be seen from the figure 2 and figure 3 that the bounds on Mx/gx necessarily 
rules out most of the parameter space in gx — M x or gx — u plane which give correct dark 
matter properties. Only a narrow region of parameter space near the s-channel resonance 
Mx ~ 2m dm is left. If dark matter mass is light, a few tens of GeV, then additional 
bounds from LEP-II experiment will apply on neutral gauge boson and its coupling. The 
agreement between LEP-II measurements and the standard model predictions forces the 
mass of additional neutral boson to be greater than 209 GeV or the couplings to be smaller 
than or of order ICE 2 [22], This will further reduce the parameter space to the region with 
gx < io~ 2 . 

After Ending the parameter space which keeps the total abundance of Nr and S 2 R within 
the Planck limit on dark matter abundance (1), we also show the relative contribution of Nr 
and S- 2 R to dark matter relic abundance in figure 4. It can be seen from the figure that Nr 
can give rise to 26 — 28% of dark matter relic density whereas S 2 R gives rise to the rest of it. 
Since their relative abundances are different, their scattering probability at direct detection 
experiments will also be different. We have taken that relative factor into account while 
calculating the dark matter direct detection cross section. 

In order to fit our model with the observed 3.55 keV X-Ray line data [1], we follow the 
constraint on the decay width of the heavier dark matter candidate Nr as obtained in[8] 

r (Nr ->■ 5 2r7 ) » 6.2 x IQ- 47 rn Nlt GeV, (18) 
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FIG. 5: Radiative decay of Nr into S 2 r and a photon 7 . 

where Nr contributes around 50% to dark matter relic abundance. Here the dependence on 
arises via the number density of dark matter. In our model, the heavier dark matter 
particle Nr can decay into the lighter dark matter particle S 2 r and a photon 7 only at two- 
loop level. The corresponding Feynman diagrams can be seen in figure 5. In the following we 
try to make a rough estimate of the decay width due to these two-loop feynmann diagrams, 
assuming all the fields inside the loop to be around TeV scale. Using our estimate for this 
decay width expressed in terms of the couplings and also using above constraint (18) we 
obtain conservative bound on the product of various couplings as discussed below. 

As Nr and S 2 r are Majorana neutrinos, corresponding to all diagrams in figure 5 there 
will be conjugate diagrams where photon connects to opposite sign particles with respect to 
diagrams in figure 5. Considering those conjugate diagrams and also considering all heavy 
masses in the internal line of the diagrams almost degenerate and neglecting electron mass, 
the decay width can be approximated as 


T(Nr S 2R 7 ) ~ 


™n. 


R ~ m S 2R 


167T m% R 


m 2 NR -m 2 S2R ) 2 [F 2 + F 2 ] 


(19) 


where 


F 1 ~ 2 


m N „ - rn s , 2 


rru- < 


I ; F 2 ~ 2 


m N „ + m S2 




I ;/ ~ 


hNfef 12 V u 5 
2567r 4 ?n^- 


( 20 ) 


Here, F\ and F 2 correspond to Lorentz invarient form factors connected to electric dipole 
moment transition and purely magnetic moment transition respectively. In general both of 


16 












them will be contributing to such decays if there is CP violating interactions involving Nr as 
well as S 2 r. However, here we assume that CP is not violated by the interactions involving 
Nr and S 2 r. Then two possible CP eigenvalues G {+i, — i} are possible for Majorana 
particles like Nr and S 2 r. So, either Nr and S 2 r will have same CP eigenvalues or opposite 
CP eigenvalues. For same CP eigenvalues we get purely electric dipole moment transition 
for which F 2 = 0 in eq. (20) and for opposite CP eigenvalues we get purely magnetic 
moment transition [24, 25] for which F\ = 0. Taking the constraint (18) into account, we 
have 

~ 6.2 x 10 ~ 47 M Nr GeV (21) 

where 

A k = 7 jin r — ms 2R = 3.55 x 10~ 6 (GeV). (22) 

Taking m NR « ms 2R ~ 35 GeV, m^- ~ 500 GeV and u 5 ~ 5 TeV, this constraint can be 
naturally satisfied if 


^Nf&f 12 V ~ 1 (for same CP eigenvalues). 

~ 10 -14 (for opposite CP eigenvalues). (23) 

As will be discussed in the next sections, to satisfy the constraints from light neutrino masses 
as well as the one-loop mass splitting between N r and S 2 r, the opposite CP eigenvalues of 
Nr and S 2 r seem to be appropriate. 

After constraining the model parameters from dark matter, collider as well as the 3.55 
keV X-ray line data, we check if the model can explain the galactic center gamma ray excess 
for the same region of parameter space. Recent analysis [15] of the Fermi Gamma Ray 
Space Telescope data has shown an excess of gamma rays with a peak of 1 — 3 GeV in the 
region surrounding the galactic center. Also reported by earlier analysis [26], the spectral 
shape of the gamma rays has a feature which resembles annihilating dark matter. One 
obtains better fit of the observed gamma ray excess from the annihilation of dark matter to 
bb pairs for dark matter mass 35 GeV[15] with the constraint on the cross-section ( av) = 
(0.77 — 3.23) x 10 -26 cm 3 /s[20]. Fit with the observed gamma ray excess for other masses 
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of the dark matter with different annihilation channels has been discussed in reference [15]. 
In our subsequent analysis we have considered particularly dark matter mass 35 GeV. Since 
we have two dark matter candidates with a mass difference of 3.55 keV, we consider the 
mass of Nr to be 35 + 3.55 x 10 -6 GeV and that of S 2 R to be 35 GeV. We then show 
in figure 2 and 3 the region of parameter space in gx — M\ and gx — u planes for which 
the total annihilation cross section of Nr, S 2 R matches the one mentioned above in order 
to produce the observed gamma ray excess. We include their relative abundance factors 
while calculating the annihilation cross sections needed to produce galactic center gamma 
ray excess. The allowed mass of neutral vector boson becomes around 70 GeV which faces 
severe constraints from LEP-II data and further constrains the coupling gx < 10 2 , as 
discussed earlier. Similarly, we find small allowed parameter space for gx < 10~ 2 and u > 2 
TeV where for simplicity we have assumed u\ = U2 = U3 = U4 = u as mentioned just above 
equation (13). 

It should be noted, taking into account of the recent work [27], that the model of one- 
loop radiative neutrino mass with dark matter originally proposed by Ma and popularly 
known as "scotogenic" model [13] suffers from a hierarchy type problem which can spoil 
the phenomenological success of this model. The problem occurs due to the contributions 
from heavy Majorana neutrinos through renormalisation group evolution (RGE) to the mass 
parameters of the scalar fields which are odd under the unbroken Z 2 symmetry of the model. 
If, for some parameter space of the model, the mass parameters of the Z 2 - odd scalars turn 
negative at some energy scale, it will break the Z 2 symmetry resulting in the loss of a cold 
dark matter candidate in terms of the lightest Z 2 -odd particle. Although the model we are 
studying is an extension of the original Ma’s model by a gauge symmetry U(l)x, the effective 
low energy model below the U{l)x breaking scale is a Ma type model with heavy Majorana 
neutrinos and a Z 2 symmetry under which several scalar fields are odd. As shown by the 
authors of [27], one can prevent the mass parameters of the Z 2 - odd scalars from turning 
negative if the physical Z 2 -odd scalar masses are restricted to a certain ranges, typically of 
the order of the heaviest Majorana neutrino. Since our dark matter candidates are singlet 
Z 2 -odd Majorana fermions and instead of Z 2 -odd scalars, these constraints can be satisfied 
easily without affecting the dark matter phenomenology discussed above. 
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IV. LIGHT NEUTRINO MASS 


The origin of tiny neutrino masses was discussed in details in [11], The tiny masses can 
arise at both tree level as well as one-loop level through the Feynman diagram shown in 
earlier works [11, 12, 14], Since out of the three singlet neutrinos Nr, Sir, S 2 r, only Sir 
gives rise to a Dirac mass term mo = yv 1 for the neutrinos through the vev of $1 (denoted 
as Vi), only one of the neutrinos acquire a non-zero mass at tree level through type I seesaw 
mechanism [28]. The tree level mass for the light neutrino in terms of the Dirac mass term 
and the mass of the heavy singlet neutrino Sir ( Mg 1R = fsUi ) can be written as 

2 y 2 v\ 


rn,. 


fsu 1 


(24) 


From figure 3, we see that the allowed region from dark matter as well as collider constraints 
suggest U\ = u 2 = U 4 = u 5 = u > 2 TeV, so we have taken u\ ~ 5 TeV and f s u\ ~ 2 TeV. 
Since V\ ~ 100 GeV, for light neutrino masses to be at sub-eV scale, the equation (24) 
suggest that the Yukawa couplings y have to be around 3 x 10 -6 which is approximately 
same as the electron Yukawa coupling in the SM. The other two SM neutrinos can acquire 
non-zero masses only when one-loop contributions are taken into account. As discussed in 
[11], the one-loop contribution (M u )ij to neutrino mass is given by 

,f',i .f5 V\V‘2U\U4 

k 


{M v 


*? 


I 671" 2 


Y hN,E it hN,z jh (Ak + 


(25) 


Assuming all the scalar masses in the loop diagram to be almost degenerate and written as 


m sc then 


T (-Bfc)ij; l ~ FTl2k 


™ 2 sc + m lk 


m lc i m2 sc - m 


2 k) 


(2 ^£ 7 ) ^2k 1 ( 9 / 2 \ 

A -^fln {m 2 Jm 2 2k ) 


(mi 


m 


2 k) 


(26) 


where = m 2 k- For fermion singlet light dark matter, rriak m sc and hence the 

above expression can be approximated as 


Ak + (Bk) 




m 2k 

ml 


The one-loop neutrino mass can be written as 

.f:S.f5 v l v 2U\U4 

k 


(M u 


*3 


167T 2 


h N^ik h NX ]k 


m 2 k 
ml 


(27) 
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Taking iq, w 4 to be at 5 TeV and m sc ~ 500 GeV, v\, v 2 at electroweak scale and the singlet 
mass rri 2 k at 100 GeV, the above expression can give rise to eV scale neutrino mass if 

/3/5^/v^iv ~ 4 x 10 11 

whereas for singlet mass m 2 k at 35 GeV, this constraint becomes 

hhh N h N ~ 1.12 x 10- 10 (28) 

One may note here that the appropriate explanation for the Gamma ray excess from the 
galactic center can be explained if one considers mx R ~ ms 2 R ~ m 2 k — 35 GeV [20]. 

The expressions for neutrino masses and mixing angles can be derived using the param¬ 
eters of the model which appear in the Lagrangian written for the model with additional 
U(l)x gauge symmetry. The value of these parameters change while going from high energy 
scale with unbroken U(l)x down to electroweak scale due to the effects of RGE. The cor¬ 
responding changes in the neutrino parameters under RGE for the original Ma model were 
studied by [29]. The authors have shown that the RGE effects on neutrino parameters can 
be quite large in these models, due to the dependence of neutrino parameters on products 
of several couplings of the model when neutrinos obtain masses only via one-loop diagrams. 
However, in our case, there is tree level contribution also to the neutrino mass matrix and in 
renormalising the theory the counter-terms of the broken electroweak phase can be obtained 
from the symmetric phase by simple algebraic relations [30]. So the analysis done in [29] is 
not directly applicable in our case. However, a more complete analysis of neutrino masses 
and mixing in our model should include RGE effects. 

V. MASSES OF DARK MATTER N R AND S 2 r 

For the discussion on dark matter, we assumed the mass of the heavier dark matter 
candidate N R to be 35 GeV, while keeping the lighter dark matter mass lower by about 3.55 
keV. From the Yukawa Lagrangian of the model in equation (3), it can be seen that N R and 
S 2R receive tree level masses from the vev’s of Xi and y 2 respectively. Considering Yukawa 
couplings fs 2 and fx of the order of 1.7 x 1CT 2 and the vev’s of y 2 and Xa of the order of 2 
TeV, with slight difference in the values of these Yukawa couplings, it is possible to generate 
mass differences of 3.55 keV between N R and S 2R at the tree level. 
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However such small mass difference at the tree level could change if there is corrections 
from higher order or loop effects. Here we have checked whether the higher order one-loop 
diagram as shown in figure 6 could potentially contribute and change the tree level estimate 
of the mass difference between S 211 and Nr. The mass matrix Mdark in the Nr, S- 2 r basis 
can be written as 


Mdark 


/aTU4 Mdarkl2 
MdarkZi fs2U2 

The one-loop contribution Mdark 12 = Mdark* 2 i to 2 x 2 mass matrix is given by 


(29) 


Mdark 12 


^ 2y 2 y*h N f 12 f 5 *i>lv2U4 

I 67 r 2 

I ’ ^X3 R1 I ’ ^X3 1 J 


(30) 


in which 


I (a, a, b, c ) 


a 2 ln(a 2 /c 2 ) — a 2 + c 2 


a 2 (a 2 — c 2 ) 


2 ’ 


(31) 


for b « a, c and 

I (a, b, c, d ) 


a 2 — b 2 


a 2 -d 2 




(32) 


for c << a,b,d and these limits are useful as active neutrino mass scale m v is very small 
in comparison to other masses. Here, m^o and are the masses corresponding to 
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Re[(p 3 ] and Im[(j) 3 ] respectively whereas m X3R and m X3I are the masses corresponding to 
Re [x’ 3 ] and Im[xl\ respectively. As denoted in the previous section, v % and iq are the vevs’ 
of electroweak scalar doublets and additional singlet scalar fields respectively. We have 
considered m v ~ and m s 1R fsU 1 - 

Particularly, if we consider m^o ps m X3 and denote the near-equal masses of Re[m s 3 ] and 
Im[m s s] as rn s where S 3 = 03 , X 3 , then 


M, 


dark 12 


2 y' 2 y*h N f 12 f 5 * v\v 2 u A 


I 67 r 2 


<t>3R <t>31 


{In K/(/M) “ 1 } 


- fs u i 




(33) 


From tree level neutrino mass ~ ^ , if we keep the mass of singlet neutrino mg 1R ~ / S W! 

hxed at 2 TeV or so as discussed just after eq. (24), then the Yukawa couplings have to 
be of the order of 3 x 10 ~ 6 to give rise to neutrino mass of order 0.1 eV. Considering 
m s = (m ^+ m ^)/2 ps 500 GeV, « 10 GeV and other parameters mentioned 

above we obtain 


M dark i 2 ~ hfi2h N x 10~ 16 GeV (34) 

Even if we consider the couplings / 5 ,/i 2 and /i,y of 0(1) then also M darkl2 ~ 10 -16 GeV 
which is O(10 10 ) less than the mass splitting of 3.55 keV as considered at the tree level. 
Therefore, with our choices of parameters, M darkl2 is negligible in (29) and hence the higher 
order one-loop correction does not alter the tree level mass difference of Nr and S 2 r- 
To check the compatibility of different requirments associated with 

1 . the condition on the product of couplings from the decay of the heavier dark matter 
producing the 3.55 keV line. 

2 . obtaining appropriate heavier neutrino mass from the tree level. 

3. obtaining appropriate lighter neutrino mass from tree level as well as one-loop level. 

4. keV mass difference of the two dark matter Nr and S 2 r at the tree level. 

we analyse the conditions in equations (23), (24), (28) and (34) respectively. Although the 
equation (34) is trivially satisfied with any value of Yukawa couplings less than 0(1) but 
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other equations give certain conditions on the Yukawa couplings as well as conditions of 
CP eigenvalues of the two dark matter particles. From equation (24), with our appropriate 
choice of parameters we have obtained y ~ 3 x 10 -6 . Considering this value of y we find that 
the condition given in equation (23) can be satisfied only for the opposite CP eigenvalues 
of two dark matter particles as other Yukawa couplings in that equation are expected to 
be of 0( 1) or less. Considering the condition of opposite CP eigenvalues in equation (23) 
and using the value of y as obtained from equation (24) we get ~ jjpZ using which in 
equation (28) we obtain the following condition on the ratio of different Yukawa couplings 

pf 2 

Mil _ 10 —6 ( 35 ) 

J3/5 

Considering Yukawa couplings like / 6 , /i 2 , and / 3 of the same order they turn out to 
be 0( ICC 3 ) whereas ~ 0( 10~ 2 ). So, in our model with suitable values of the Yukawa 
couplings as discussed above, it is possible to explain the 355 keV X-ray line, tiny neutrino 
mass and gamma ray excess from galactic centre in a single framework. 

VI. RESULTS AND CONCLUSION 

We have studied an abelian gauge extension of the standard model which can predict tiny 
neutrino mass and stable dark matter candidate naturally. In particular, we have discussed 
the possibility of explaining the recently observed 3.55 keV X-ray line and the galactic center 
gamma ray excess from a common dark matter origin within the framework of this abelian 
gauge model. Although the model has both scalar and fermionic dark matter candidates, we 
choose to study only fermionic dark matter candidates to serve our goal better. The dark 
matter candidate is guaranteed to be stable by a remnant Z 2 symmetry after the abelian 
gauge symmetry gets spontaneously broken. In order to explain the 3.55 keV X-ray line, we 
assume the dark sector to consist of two dark matter particles: the lightest Z 2 -odd particle 
(S 2 r) and the next-to-lightest Z 2 -odd particle Nr, both of which are singlet fermions. The 
mass difference between the two dark matter particles is chosen to be 3.55 keV such that the 
heavier one can decay into the lighter one and a photon at loop level. In order to explain 
the galactic center gamma ray excess, we choose the lightest dark matter mass to be 35 
GeV and check whether the two dark matter particles give rise to the required annihilation 
cross sections. We also take into account the constraints from dark matter direct detection 
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experiments like XENONIOO on spin dependent scattering cross section of dark matter off 
nuclei. These models can also face stringent limits on new gauge boson mass M x and gauge 
coupling gx- Using the results from [23] where the authors found the lower bound on M x /g x 
to be 2.6 TeV for BR(X —>■ SM) = 90% and gx — 0-1 we also use three different cuts on 
Mx/gx starting from a moderate 2 TeV to a conservative 4 TeV on Mx/gx■ These limits 
will be even weaker in those region of parameter space where gx can be much lower than 
0.1. We find that, even after applying a conservative lower limit on M x jg x as 4 TeV, we 
still have some parameter space available near the s-wave resonance region which can satisfy 
all constraints related to dark matter and colliders. 

After showing the allowed parameter space in terms of gx,Mx as well as u , the com¬ 
mon vev of the scalar singlets, we constrain the other parameters of the model from the 
requirement of producing the correct 3.55 keV X-ray flux, sub-eV neutrino mass and keeping 
one-loop mass splitting between two dark matter candidates below keV. From the X-ray 
flux constraints, we find that the product of four relevant dimensionless couplings have to 
be around 1 or 10“ 14 for same CP eigenvalues or opposite CP eigenvalues respectively for 
Nr and S 2 r- Similarly, the constraints from sub-eV neutrino masses keep the product of 
four relevant dimensionless couplings tuned at around 10~ 10 for singlet fermion dark matter 
masses of a few tens of GeV. Constraints from tree level light neutrino mass and decay of 
heavier dark matter to lighter dark matter and photon show that only opposite CP eigen¬ 
values of N r and S 2 r could be possible for Yukawa couplings of 0( 1) or less. In our model 
the small mass difference between two dark matters of 3.55 keV considered at the tree level 
remains unchanged even with higher order corrections. Also the allowed parameter space 
in gx — M\ plane is very limited. This is because for light dark matter mass, in order to 
explain GC excess and 3.55 keV X-Ray line together, the constraints from dark matter exper¬ 
iments as well as bound on M x /gx allow only a limited region near the s-channel resonance 
M x ~ 2771dm- For dark matter mass around 35 GeV, the allowed mass of neutral boson 
becomes around 70 GeV which again faces severe constraints from LEP-1I data and further 
constrain the coupling g x < 10” 2 . Due to the very limited parameter space available, this 
model will undergo serious scrutiny at future experiments with more sensitivity. 
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